Astronomy & Astrophysics manuscript no. radio'3 


February 2, 2008 


(DOI: will be inserted by hand later) 





Millimeter-wave emission during the 2003 low excitation phase 

of rj Carinae 

Z. Abraham^, D. Falceta-Gongalves^ , T. P. Dominici^, L.-A. Nyman^''^, P. Durouchoux^, F. McAulifFe^, A. 

Caproni^, and V. Jatenco-Pereira^ 



in 
o 
o 

(N 

Ph. 
< 



> 
o 

00 

O 
in 
o 

:^ 

6 



X 



^ Institute de Astronomia, Geofi'sica e Ciencias Atmosfericas, Universidade de Sao Paulo, Rua do Matao 1226, 

Cidade Universitaria, 05508-900, Sao Paulo, Brazil 

e-mail: zulema@astro.iag.usp.br 
^ European Southern Observatory, Casilla 19001, Santiago 19, Chile 

Onsala Space Observatory, 439 92 Onsala, Sweden 
* Service d'Astrophysique (CNRS URA 2052), CEA Saclay, 9119Gif-sur-Yvette, France 



Received; Accepted 

Abstract. In this paper we present observations of rj Carinae in the 1.3 mm and 7 mm radio continuum, during the 
2003.5 low excitation phase. The expected minimum in the light curves was confirmed at both wavelengths and 
was probably due to a decrease in the number of UV photons available to ionize the gas surrounding the binary 
system. At 7 mm a very well defined peak was superimposed on the declining flux density. It presented maximum 
amplitude in 29 June 2003 and lasted for about 10 days. We show that its origin can be free-free emission from 
the gas at the shock formed by wind-wind collision, which is also responsible for the observed X-ray emission. 
Even though the shock strength is strongly enhanced as the two stars in the binary system approach each other, 
during periastron passage the X-ray emission is strongly absorbed and the 7 mm observations represent the only 
direct evidence of this event. 
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1. Introduction 

Even after more than a century of intensive studies, lit- 
tle is known about Eta Carinae as a star, hidden as it is 
by a dense cloud of gas and dust. Analysis of the histor- 
ical optical data revealed a 5.52 year periodicity, charac- 
terized by very sharp dips in the light curve, lasting for 
only a few months (Damincli 1996). Similar periodicity 
was found in the radial velocity of the Pa7 and Pa,S lines, 
suggesting the existence of a binary system in a highly 
eccentric orbit (Damineli, Conti & Lopes 1997; Davidson 
1997; Damineli et al. 2000), although the orbital param- 
eters were not well defined by the observations. Damineli 
et al. (2000) derived an eccentricity 0.65 < e < 0.85 from 
the radial velocity data and Corcoran et al. (2001a) found 
an even higher value, e = 0.9, in a model where the X- 
ray emission is produced by wind- wind collisions. The 
ephemerides for periastron passage was also constrained 
by the X-ray observations, since the strong periodic in- 
crease in luminosity was attributed to an increase in the 
shock strength as the two stars approached each other. 
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Unfortunately, the maximum in the X-ray emission is not 
observable, because at periastron the shocked region seems 
to be absorbed by a large H column density, which pro- 
duces the pronounced dip in the X-ray light curve and puts 
constr ains on the ephem erides f or conjunction and oppo - 
sition ijCorcoran et al. 2001a: Pittard fc Corcoran 2002() . 
The existence of a binary system was contested by 
Davidson et al. (2000), based on high spatial resolu- 
tion HST spectroscopic observations. Recently Falceta- 
Gongalves, Jatenco-Pereira & Abraham (2005) gave an al- 
ternative explanation to the X-ray absorption, which puts 
the periastron close to conjunction, in agreement with the 
HST data. 

Radio observations of 77 Carinae, from 6 cm to 
1.3 mm, show also the periodic, sharp dips in their 
light curves (Cox et al. 1995a; Duncan, White & 
Lim 1997; Abraham & Damincli 1999). The continuum 
emission is attributed to thermal bremsstrahlung from 
an optically t hick plasma, pr oduced either in the rj 
Carinae wind lICox et al. 1995al) or in an equatorial disk 
l|Duncan k White 20031)" Interferometric monitoring with 
the Australian Telescope Compact Array at 3 and 6 cm 
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Fig. 1. mm- wave hght curves of r] Carinae at the epoch 
of the predicted low excitation phase. Open circles corre- 
spond the 1.3 mm emission (right axis), full circles corre- 
spond to 7 mm emission (left axis). 

with 1" resolution, showed an elongated structure, reach- 
ing a maximum size of about 4" (1.3x 10^^ cm at a distance 
of 2.3 kpc), which shrink s to less than 1" dur ing the min- 
imum in the light curve ijPuncan et al. 1997j) . This effect 
was interpreted by Duncan & White (2003) as due to a 
decrease in the number of ionizing photons during the low 
excitation phase. 

In this paper we report contemporaneous observations 
of the last minimum in the light curves, which started 
around May 2003, at two radio wavelengths: 1.3 mm and 
7 mm, obtained respectively with SEST, at La Silla, Chile 
and with the Itapetinga radiotelescope, in Brazil. The 
closely spaced observations provided not only the light 
curves with great precision, but also what we believe was 
the detection of radio emission from the colliding winds, 
highly enhanced during periastron passage. Since at that 
time the X-ray emission from the winds was totally ab- 
sorbed, our observations represent the first detection of 
this phenomenon. 

2. Observations and Results 

To better understand the radio light curve, we made daily 
observations at 7 mm with the Brazilian Itapetinga ra- 
diotelescope and weekly observations at 1.3 mm with the 
ESQ (European Southern Observatory) SEST radiotele- 
scope at La Silla, Chile, during the periodic event pre- 
dicted for June 2003. 

The Itapetinga radiotelescope has a beam size of about 
2 arc-min at 7 mm, enough to separate rj Carinae from 
the other HII regions in the Carll cloud; the observations 
consisted of scans across the source with 20 arc-min ampli- 
tude, to allow the subtraction of the underlying sky con- 
tribution. The atmospheric opacity was calculated using 



the technique developed by Abraham & Kokubun (1992); 
the flux density was calibrated by observation of the stan- 
dard point source Virgo A; the stability of the system was 
continuously checked by intercalating observations of the 
HII regions G287.57-0.59 and G287.55-0. 63, part of the 
extended Carina Nebulae l|Retalack 1983(l . with those of 
T] Carinae. 

Observations at 1.3 mm were made with SIMBA 
(SEST Imaging Bolometer Array), a 37-pixel bolometer 
array with a beamsize of 23" and separation between pix- 
els in the sky of 44" . The observations were made in fast 
scanning mode (80" s""'^ in elevation). The peak flux of rj 
Carinae was determined by fitting a two-dimensional gaus- 
sian to the source; series of 4-8 maps were obtained on each 
occasion, the uncertainty in peak flux was obtained from 
the rms value of the fluxes in each map. The atmospheric 
opacity was determined from sky dips and daily observa- 
tions of Uranus provided absolute calibrations. The last 
observation coincided with the end of the operations of 
the SEST radiotelescope. 

The data, presented in Figure 1, confirms the appear- 
ance of the predicted dip, with a decrease in flux density 
by a factor of four at both wavelengths. The closely spaced 
observations revealed also a short duration peak at 7 mm, 
reaching the maximum flux density in June 29 and last- 
ing for about 10 days, superimposed to the fast decline 
in emission. The peak is not conspicuous at 1.3 mm, al- 
though the time sampling at this wavelength is not as 
dense as at 7 mm. Although the general behavior of the 
light curves at radio frequencies is similar to what was 
observed in the previous cycle, made with coarse tem- 
poral resolution (Abraham & Damineli 1999; Cox 199^), 
the short-lived peak observed at 7 mm was unexpected 
and difficult to explain in terms of, for example, a sudden 
increase in the number of ionizing photons. We thus inves- 
tigated if the excess radio emission could be produced in 
the thin shock region that is responsible for the observed 
X-ray emission. 



3. Model for the mm-wave emission 

To interpret quantitatively the mm-wave light curves 
we will assume that the free-free emission is pro- 
duced both in an e xtended disk surrounding ij Carinae 
( Duncan et al. 199^ and in the coofing shock region 
formed by wind-wind collision. We will further assume 
that the flux of UV radiation that maintains the disk ion- 
ized, suddenly decreases, or maybe even stops, as proposed 
by Duncan & White (2003). 

3.1. Emission from tlie optically thick disk 

If the ionizing flux decreases, part of the ionized material 
will recombine at a rate r^cc given by: 

n-oc = a'^^ric, (1) 
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where rig is the electron number density and a^^-* the re- 
combination coefficient, given by: 



v(l) 



2.06 X 10-"Z2 



*(T/Z2) cm^ s" 



(2) 



with the function ^'(T/Z^) tabulated in Spitzer (1978). In 
the inner region of the disk, electron densities Ue ^ 10^ 
cm~"^ and temperatures T ~ 8000 K are required to ex- 
plain the large flux densities and optically thick contin- 
uum radio spectrum, as well as the HSOa, H35q; and 
H40q; reco mbination lines, which present strong ma ser am- 
plification l|Cox et al. 1995bl:lAbraham et al. 2002|) . Using 
these values in equations (1) and (2) we find a recom- 
bination timescale of 2.5 days. Since the disk density 
decreases with dista nce to t] Carinae ijCox et al. 1995at 
lAbraham et al. 20021) . the cooling time will be larger in 
the outer parts of the disk, explaining qualitatively the 
duration and shape of the declining part of the light curve. 
The net result of this process is the formation of a hole 
of neutral material in the center of the ionized disk. Since 
the actual electron density distribution in the disk is not 
known and we are mainly interested in the residual emis- 
sion, we adjusted the 7 mm flux density emitted by the 
disk by a smooth function of the ephemerides t: 



'S'disk(7 mm) = a cxp {bt + c) + mt + n. 



(3) 



The parameters a = 3.1 x 10"^ b = -0.228, c = 7.50, 
m = 0.004 and n — 0.94, were determined from the data 
at both sides of the peak, the time t was measured from 
29 June 2003 (JD = 2452829). 

In Figure 2 we show an expanded version of the 7 mm 
light curve, with the assumed disk flux density indicated 
by a continuous line (upper panel). The residual emission, 
which origin we will try to establish, is seen in the lower 
panel. 

3.2. Free-free emission from the ionized winds 

Let us assume that the winds are completely ionized and 
their temperatures are close to 10^ K; their densities p at 
a distance r from the corresponding star are related to 
their mass loss rates M and terminal wind velocities v by: 



M 



(4) 



The flux density at mm wavelengths must be calcu- 
lated taking into account radiation transfer since, at the 
densities and temperatures found close to the stars, the 
optical depth can be large. Using the Rayleigh-Jeans ap- 
proximation we have: 



S{X) 



2kT 



(l-e-^(^))f]. 



(5) 



where is the solid angle subtended by the emitting re- 
gion and t(A) the optical depth. 

To explain the observed flux density at the peak of the 
emission light curve (about 1 Jy), in the most favorable 
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Fig. 2. Expanded view of the 7 mm light curve. Upper 
panel: observed 7 mm flux density (filled circles) and 
model disk emission (continuous line). Lower pannel: 
residuals of the observed flux density after subtracting 
the disk contribution. Zero in the time axis corresponds 
to June 29, 2003 (JD = 2452819). 



condition of r ^ 1, it is necessary a source subtending a 
solid angle of 6 x 10~^^ srd, which at the distance of 77 
Carinae corresponds to a radius ro = 5.3 x 10^^ cm. 

The optical depth of the emitting region is found from 
the integral: 



t(A) = / K//(A)ds, 



(6) 



where Kff{X) is the free- free absorption coefficient, which 
at radio wavelengths can be calculated from: 



KffiX) = 2 X 10 



_23 n,neg//(A)A^ 

r3/2 



(7) 
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where 5// (A) is the Gaunt factor for free-free emission. 
For the physical conditions in the winds, T < 10^ K and 
hc/X <C kT, gff{X) can be approximated by: 



9ffW 



V3 



17.7 + hi( 2.2^^ 



(8) 



For T > 10^ K: 



,,,(A)^-ln 2.2— . 



(9) 



Taking Arg as the depth of the shell that gives r = 1, 
wc find from equations (4), (6), (7) and (8) Ar = 3.7x 10^^ 
cm. 

If we assume that the peak of the 7 mm emission is 
due to an increase in the mass loss rate of the star, the 
perturbation will reach the radius r — 5 x 10^^ cm in 
approximately 2 years. Since the observed time scales are 
much smaller, we can eliminate the winds as the emission 
source. We investigate now the possibility of the emission 
being produced by the wind-wind collision. 
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Fig. 3. Time evolution of the temperature and density 
behind the shock, due to energy losses by free-free and 
line emission, assuming pressure equilibrium. The initial 
conditions are: r(0) = 2 x 10** K and n(0) = 6x10^ cm'^ 



3.3. Free-free emission from shocked material 

The physical conditions in the shock produced by wind 
collision were calculated analytically (e.g. Girard & 
Wilson 1987; Usov 1992) and applied to the rj Carinae 
binary system by Ishibashi et al. (1999) and Corcoran et 
al. (2001a). Detailed numerical simulations were obtained 
by Pittard et al. (1998) and Pittard & Corcoran (2002), 
including the effects of radiative cooling. In all cases the 
spectrum of the X-ray emission and the resulting light 
curve were satisfactorily explained, using appropriate pa- 
rameters for the orbit, mass loss rates and wind velocities. 

The radio emission from wind-wind collisions was cal- 
culate by Dougherty et al. (2003), using numerical models 
and including radiation transfer effects and different incli- 
nation angles for the orbital plane, to explain the thermal 
and non-thermal emission detected in several binary sys- 
tems (Dougherty et al. 1996; Dougherty & WiUiams 2000; 
Moran et al. 1989; Churchwell et al. 1992; WiUiams et al. 
1997; Niemela et al. 1998). 

However, close to periastron passage and because of 
the very high mass loss rate from rj Carinae (Hillier et 
al. 2001) and the high eccentricity of the binary system 
(Damineli et al. 2000; Pittad & Corcoran 2002), radiative 
cooling is very fast and the shock region becomes very nar- 
row. For that reason, the temperature, density and ioniza- 
tion profiles obtained from numerical models (e.g. Pittard 
& Corcoran 2002), in which the cell size is about 3 x 10^^ 
cm, cannot be used for our purpose. Therefore, we will 
calculate the expected radio flux density at 7 mm with 
plausible combinations of temperature, density and size 
of the emitting region. 

We first determine the physical conditions at the shock 
using the expressions presented by Usov (1992). The mo- 



mentum balance surface intercepts the line joining the two 
stars at a distance ri from rj Carinae, given by: 

D , , 

ri = r, (10) 

1 + 772 

where D is the distance between the stars and 77 — 
MsVs/ MpVp. Mp and Ms are the mass loss rates of 77 
Carinae and the companion star, Vp and Vg their respec- 
tive wind velocities. Pittard & Corcoran (2002) favored 
?7 = 0.2 in their fitting of the X-ray spectra obtained with 
Chandra. 

Shocks are formed at both sides of the contact surface. 
The density p2 and temperature T2 behind each shock can 
be derived from the values pi and Ti at the shock front: 

P2 _ (7 + l)M^ f... 
pi (7-l)M2 + 2' ^ ' 

P2T2 ^ - (7 - 1) ^^2) 

PiTi 7 + 1 

where 7 is the adiabatic exponent, M = 

_ 1 

vi {jkTi/ jirriH) ^ is the Mach number, n the molecular 
weight and mn the mass of the H atom. 

Assuming for 77 Carinae a mass loss rate of 2.5 x 10^"* 
Mq yr~^, a wind velocity of 700 km s~^ and a separa- 
tion D = 1.6 AU at periastron, co rresponding to an or- 
bit with eccentric itv e = 0.9 ( Pittard fc Corcoran 20021: 
iHillier et al. 200lh . we find a number density ni = 
pi/pma = 5 X 10-^" cm"-*. For a temperature Ti — 10^ K, 
the post-shock physical parameters will be rt2 = 2 x 10^^ 
cm^'^ and T2 = 6 x 10^ K. For a secondary star with 
mass loss rate Mg ~ 10~^ Mq yr~^ and wind veloc- 
ity Vs = 3000 km s^^, the corresponding values behind 
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the shock are 712 = 6 x 10^ cm"'^ and — 2 x 10* 
K, tem perature compatible with Chandra X-ray obser- 
vations ijSeward et al. 200itlCorcoran et al. 2001h(l . 

With these initial values, we can calculate the tem- 
perature and density time evolution from the first law 
of thermodynamics assuming, as in Mathews & Doane 
(1990), isobaric equilibrium and cooling by free-free and 
line emission: 



2,kp{t) dT{t) kT{t) dp{t) 
2fj,niH dt pruH dt 

with A is giveny by: 



= -A. 



(13) 



A = 



2.4xl0"2^r5 



61 TP 



(I + 62T9) 



niUe erg cm s ^,(14) 



where and Ue are the number density of ions and elec- 
trons, respectively, and 61 = 1.53x 10^^''', 62 = 1.25x10"^, 
p = 1.2 and q = 1.85. We did not consider heating by the 
stellar UV radiation because most of the radio emission 
is produced at temperatures larger than 10^ K, when this 
heating process is not effective. 

In Figure 3 we show the time evolution of rig and T, 
for the secondary shock and the physical conditions in 
the post-shock gas derived previously. We can see that 
the temperature decreases slowly during approximately 4 
days and then drops very fast, as the density increases and 
line emission becomes effective. For the primary shock, the 
cooling time is less than a day, since the initial tempera- 
ture is lower and the density much higher. 

As in the case of the individual wind emission, we can 
estimate the size of the optically thick region necessary to 
produce the observed peak at 7 mm. Figure 4 shows its 
radius as a function of the electron temperature, for the 
conditions found in the secondary shock. This size will 
be 4 times larger for the primary shock. We should notice 
that T cannot be much larger than one at 7 mm, otherwise 
the optical depth at 1.3 mm could also be large, in which 
case the expected flux density would be 28 Jy, much larger 
than what was observed. In the same Figure we show for 
this reason, the depth L of the emitting region that gives 
r = 1 both for 7 mm and 1.3 mm. The allowed values of 
L lie between the two lines. The dotted line represents the 
separation between the secondary star and the shock at 
periastron, which gives an upper limit to the depth and 
temperature of the material responsible for the peak emis- 
sion. The same calculation can be made for the primary 
shock, although in this case, for the same emitting area 
the emission is 16 times smaller. 

4. Conclusions 

In this paper we presented the 1.3 and 7 mm light curves 
of 77 Carinae during the 2003.5 low excitation phase. The 
expected minimum in the light curves was confirmed but 
an unexpected peak, superimposed to the decreasing fiux 
density, was observed at 7 mm. It reached its maximum 
in 29 June 2003 and lasted for about 10 days. To interpret 




Fig. 4. Radius R and depth L of the emitting region that 
reproduces the peak flux density and r = 1, respectively, 
as a function of the electron temperature. The dotted line 
represents the separation between shock and star. 

the light curves we assumed that the decreasing part of the 
emission is produced in the extended disk surrounding the 
rj Carinae binary system. We tested the possibility of the 
peak emission being produced either by the single stars 
winds or by the wind-wind shock region, also responsible 
for the X-ray emission. We calculated the size of the emit- 
ting region assuming single wind emission and obtained a 
radius of 5 x 10^^ cm and a depth about 3.7 x 10^^ cm 
for T — 1. Considering the wind velocity of 700 km s^^, 
a change in mass loss rate would take about 20 years to 
reach this distance. An alternative was to investigate the 
emission from the wind-wind collision. We calculated the 
physical conditions in the material behind the primary 
and secondary shocks and estimated the cooling time due 
to free-free and line emission. The size of the emitting re- 
gion was determined as a function of temperature, as well 
as the possible depth range compatible with both 7 and 
1.3 mm observations. 
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